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Abstract 

The  ultrasonic  absorption  of  the  five  alkali  metal  sulfates  was 
measured  in  water  in  the  frequency  range  from  25-520  MHz.  Only  one  relaxa¬ 
tion  was  found  for  all  the  salts  and  in  all  systems  it  was  in  the  region  of 
100-400  MHz.  The  amplitude  of  the  absorption  increased  with  increasing  con¬ 
centration  while  the  relaxation  frequency  decreased  or  remained  constant. 

The  amplitude  and  relaxation  frequency  behaved  in  a  similar  manner  when  the 
temperature  was  decreased.  These  results  were  analysed  in  terms  of  a  three 
state  mechanism  of  association  where  the  two  reaction  steps  are  coupled. 

The  experimental  results  included  a  B'  parameter  which  was  the  residual 
absorption  due  to  a  higher  frequency  relaxation  and  solute -solvent  inter¬ 
actions.  The  three-state  mechanism  yielded  forward  rate  constants  at  0.5  M 
and  25°  C  for  the  formation  of  the  inner-sphere  complex  of  1.0  x  109  sec  1 
for  all  salts  except  potassium  sulfate  which  was  2.0  x  109  sec  1  This 
mechanism  successfully  explained  the  decrease  in  the  relaxation  frequency  of 
sodium  sulfate  with  increasing  concentration.  The  volume  change  for  the 
reactions  can  be  calculated  and  they  exhibit  a  temperature  and  concentration 
dependence  similar  to  the  absorption  amplitude.  A.t  the  higher  concentrations 
it  was  possible  to  treat  the  two  steps  as  decoupled  to  a  good  approximation . 
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Intrcduct ion 

The  alkali  metal  ions  are  extremely  important  constituents  of  the  sea1 
and  of  biofluids. p  In  the  past  the  characteristic  kinetics  of  reaction  of 
these  ions  in  aqueous  solution  were  considered  immeasurably  fast.  With  the 
advent  of  relaxation  methods  Eigen  and  Maass3’4  were  able  to  investigate  the 
kinetics  of  association  of  these  simple  metal  ions  with  the  complex  multi- 
dentate  ligands  NTA3”,  EDTA4-  and  DSLTA4".  In  addition  other  kinetic 
studies4’5  have  been  carried  out  on  the  interaction  of  the  alkali  metal  ions 
with  other  large  molecules,  particularly  in  methanol  where  association  is 
very  extensive.  However,  no  systematic  studies  have  been  carried  out  on  the 
kinetics  of  association  of  alkali  metal  ions  with  simple  anions.  We  were 
particularly  interested  in  the  association  with  sulfate  because  of  the 
importance  of  sulfate  association  in  seawater  and  the  extensive  studies  made 
of  sulfate  association  with  +2  and  +3  cations.6’7’8’9 

The  problems  with  the  study  of  alkali  metal  ion  -  sulfate  kinetics  in 
water  are  two-fold.  First,  the  association  constants  are  very  small11  as 
are  the  important  thermodynamic  parameters  AH  and  AV.  Secondly,  the  rates 
are  expected  to  be  very  fast.10  For  these  reasons  we  decided  to  use  ultra¬ 
sonic  absorption.  Preliminary  results  showed  measurable  effects  although 
it  was  clear  very  early  that  we  could  expect  an  unpleasant  combination  of 
very  short  relaxation  times  and  very  low  absorption  amplitudes. 

Eigen  and  Maass3  used  the  following  mechanism  to  treat  the  results 
of  their  alkali  metal  work: 

1  k  11  k?  111 

M+(aq)  +  Xn_(aq)  *  [m+~ -Xn-]^ aq^  =  M+Xn“(aa)  (l) 

k  k 

“1  -2 


- 1 
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State  I  is  the  separated  aqueous  ions;  state  II  an  outer-sphere 
complex;  and  state  III  nn  inner-sphere  complex.  The  fact,  that  only  one 
relaxation  was  observed  for  these  systems  makes  the  above  the  most  appropriate 
scheme.  The  characteristic  water  substitution  rates3  for  the  alkali  metal 
ions  are  all  on  the  order  of  10°  sec  1  so  slow  relaxations  are  not  to  be 
expected  in  these  systems.  The  four  state  model  of  Eigen,  Diebler  and  Tamm3  ‘ 
seems  overly  complex  for  the  alkali  metal  ion  reactions.  In  this  model  step 
I  is  a  diffusion-controlled  approach  of  the  aquo  ions  while  steps  II  and  III 
represent  the  step-wise  loss  of  waters  from  anion  and  cation  to  give  the 
contact  ion  pair.  For  metal  ions  whose  characteristic  water  exchange  rates 
are  fairly  slow  the  three  steps  can  be  observed.  For  alkali  metal  ions  with 
their  high  water  exchange  rates,  the  steps  are  too  close  in  rate  to  be 
distinguishable.  Thus,  the  three-state  model  seems  to  be  most  appropriate 
to  our  case.  The  first  step  is  the  diffusion  controlled  approach  of  the 
aquo  ions  and  the  second  is  the  slower  formation  of  the  inner-sphere  complex. 

Inner-sphere  complexes  of  LiS04  and  NaS04  have  been  observed  by  IR 
spectrophotometry  in  fairly  concentrated  solutions.13’  14  Other  Raman 
studies31’  have  indicated  that  sodium  sulfate  forms  only  an  outer-sphere 
complex.  However,  the  Raman  studies  were  done  at  lower  concentrations  than 
either  the  IR  or  our  ultrasonic  studies.  Therefore  the  concentration  of 
inner -sphere  complexes  was  probably  too  small  to  detect. 

Experimental 

Materials:  Lithium  sulfate  monohydrate,  anhydrous  sodium  sulfate  and 
anhydrous  potassium  sulfate  were  Fisher  Certified  ACS  reagents  and  were 
used  as  received.  Cesium  sulfate  was  Fisher  reagent  grade  and  was  used  as 
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received.  Rubidium  sulfate  •..•as  prepared  from  silver  sulfate  and  rubidium 
chloride  and  purified.  Stock  solutions  of  all  the  suits  were  prepared 
volumetrically  using  deionised  water. 

Apparatus:  The  send-recei ve  pulse  technique  was  used  for  all  ultra¬ 
sonic  measurements.  In  this  technique  measurement  of  pulse  attenuation  as  a 
function  of  path  length  yields  the  ultrasonic  absorption  coefficient.  A 
Matec  6600  pulse  modulator  and  receiver  with  different  plug-in  units  furnished 
the  basic  electronics.  The  plug-in  units  used  were  Model  760  (10-90  MHz), 
Model  765  (90-500  Mis)  and  Model  770  ( 500-700  MHz).  The  output  signal  from 
the  appropriate  plug-in  is  gated  into  a  Matec  1255  A  pulse  amplitude  monitor. 
This  device,  basically  a  peak- reading-voltmeter,  provides  a  direct  reading 
of  the  amplitude  of  the  selected  pulse  in  decibels.  The  frequency  is 
measured  by  beating  the  pulse  with  n  Hewlett-Packard  608d  or  612A  signal 
generator,  then  measuring  the  signal  generator  frequency  at  zero  bent 
frequency  by  means  of  a  Hewlett-Packard  5527C  frequency  counter. 

Two  cells  were  employed.  The  medium  frequency  cell  uses  a  matched 
pair  of  5  MHz  X-cut  quarts  transducers  and  has  an  effective  range  of  15- 
205  MHz  in  water.  The  high  frequency  cell  employs  50  MHz  X-cut  crystals 
permanently  bonded  to  quartz  delay  rods  and  has  a  effective  range  of  90-520 
MHz  in  water.  For  the  high  frequency  cell  a  new  path-length  measuring  device 
was  developed.  We  used  a  Hewlett-Packard  displacement  transducer  that 
develops  a  voltage  directly  proportional  to  distance.  This  device  enables 
us  to  semi-automate  the  high  frequency  measurements  since  ve  car.  automatically 
process  signals  related  to  onlsc  amplitude,  path  length  and  frequency  using 
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a  computing  calcalutor.  The  medium  frequency  technique  has  been  described 
in  the  literature1*'  and  the  high  frequency  technique  will  be  submitted  for 
publication. 1_* 

The  temperature  of  the  jacketed  cells  was  controlled  to  +  0.05°  C 
with  a  Lauda  TK  j!K)-D  temperature  bath  and  monitored  with  a  Yellow  Springs 
Instrument  4^5  C  Telethennometer . 

Treatment  of  Data 


where  f  =  relaxation  frequency 

A  =  relaxation  amplitude 

B  =  high  frequency  background 

In  dilute  solutions  exhibiting  only  one  relaxation  3  is  very  close  to  the 
pure  solvent  absorption  value.  In  the  ease  of  the  concentrated  alkali  metal 
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sulfate  solutions  examined  here,  v;e  should  assume  that  the  measured 
(a/f3)  should  really  be  considered  as 


(— )  =  - - -  ♦  (  - - - — 

f2  i  ♦  (r/fr>*  V  i  +  U/r)2 


measured  observed 

relaxation 


high  frequency 
relaxation 


nnd  B  =  pure  solvent  absorption 

/SB  =  change  in  pure  solvent  absorption  caused  by  presence  of  solute. 

Then  we  shall  treat  our  data  by 

(— )  =  (— ~  -  B)  s  - - - ♦  B'  (5) 

f2  f2  1  +  (f/fr)2 

So  B'  simply  includes  any  small  amplitude  relaxation  above  our  high  frequency 
limit  plus  changes  in  the  solvent  background.  It  should  be  a  constant  for 
a  particular  salt  at  a  given  concentration  and  temperature. 

The  amplitude  of  the  relaxation,  A,  and  the  relaxation  frequency,  f  , 
are  related  to  the  relaxation  time  nnd  thermodynamic  properties  of  the  process 
occurring 
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RT*  X 


t"1  =  2n  f 
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where 


X  = 


(6) 

(T) 

(3) 


-•  stoichiometric  coefficient  of  reactant  i 

-  molar  concentration  of  reactant  i 

0  =  expansibility  of  solution 

-  velocity  of  sound  in  the  solution 
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0  -  compressibility  of  solution 

V  =  volume  of  solution 

nnd  the  other  symbols  have  their  normal  meanings  fo:1  »  chemical  equilibria'. 

The  expression  for  the  relaxation  time  in  terms  of  individual  rate 
constants  can  be  derived  from  a  consideration  o'"  the  three-state  model 
given  above.1'*  We  shall  not  assume  that  the  two  steps  are  decoupled  but 
shall  assume  that  the  activity  coefficients  are  not  changing  appreciably 
with  concentration  in  the  range  of  interest,  0.^0  to  l.Q  molar.  Since 
the  activity  coefficients  show  a  broad  minimum  with  concentration  in  the 
range,  this  latter  assumption  should  be  quite  good.1'* 

Using  sodium  sulfate  as  an  example 


II 


kj  k-? 

Na*(nq)  +  S04‘s(nq)  ~  [Ka-SO*]"  »  [n»S04]‘ 

k  k 

- 1  -2 


Then 


K£  ”  KI(1 


1  .  k  )  . 

11  [»«']  [so,-5] 


1°) 


(10) 


Kj.  =  overall  association  constant  (available  from  equilibrium  data) 
Kj  =  outer-sphere  association  constant 
Kjj  =  inner-sphere  association  constant 
[a]  =  activity  of  species  A 


We  shall  calculate  from  the  Eigen-Euoss  equation.10’*'1  can  then 

be  evaluated.  Using  the  rate  equations  for  the  formation  of  (fla •  S04 )  and 
(NaSO.. )  nnd  the  macs  balance  conditions,  v.e  can  obtain  expressions  for  the 
relaxation  tunes  using  standard  techniques. 
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The  eigenvalues  of  the  matrix  give  the  relaxation  times 
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We  will  put  this  equation  in  more  usable  form  by  using  the  experimental 

Ky>  the  theoretically  calculated  K_  anu  the  theoretical  value  for  k  *xv’ 

*“  I  “l 

In  addition  we  shall  use  Pitzer's  theoretical  equation^  for  activity 
coefficients. 

The  expression  for  isrl 
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and  for  k  is"0 
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If  the  diffusion  coefficients,  and  are  not  kr.ovn  they  can  be  estimated 
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from  the  expression 
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Robinson  and  Stokes  "’*  give  \?  as  a  function  of  temperature.  Pitrer's 
equation  for  activity  coefficients  fits  the  experimental  data  up  to  ionic 
strengths  of  5  M.  Tt  is  given  in  the  form 
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v.’here  f  ,  and  C  ^  are  explicit  functions  of  the  ionic  strength  and 

three  adjustable  parameters  that  vary  from  salt  to  salt.  The  original 
paper^5  should  be  consulted  for  further  details. 

The  combination  of  there  calculations  gives  the  following  expression 
for  the  relaxation  times 


I, II 


k  k  K '  ,  k  Ki  K ' 
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\t  v. 

since  ( — - — —)  =  V+ 

The  activity  coefficient  identity  results  from  treating  M^SO^  as  a  mixture 
of  (2I+,  S04'?)  and  (M+ ,  NtSO.,")  -..here 

Y*  (M+  )  y _  (S04”2)  =  V+3  (MpSOj 

and  Y+  (M+)  Y_  (KS0„")  =  vj  (M-S04) 

One  can  then  solve  for  V^Y_  end  v_  in  terms  of  y  which  in  then  available 
from  Pitzer's  equation.  When  K^.  is  known  all  the  rate  constants  for  the 
three  state  model  can  be  determined  once  an  'a'  value  for  the  outer-sphere 
complex  is  known.  At  the  lower  temperatures  where  values  are  not 
available,  estimated  values  must  be  used  to  obtain  (kg,  k_2)  values. 

Results  and  Discussion* 

The  excess  ultrasonic  absorption  for  different  concentrations  of 
Ka^SO-t  is  shown  in  Fig  1.  The  precision  of  the  individual  data  points  is 
very  good  for  all  solutions  measured  with  the  largest  deviation  of  any 
individual  point  from  a  derived  single  relaxation  line  being  1^.  with  the 
low  amplitude  relaxations  in  these  system^  such  precision  is  necessary  for 
meaningful  data  treatment.  As  can  be  seen  from  Fig  1  the  excess  absorption 
does  increase  with  increasing  salt  concentration.  Below  0.10  M  Na^S04 
it  is  difficult  to  detect  a  relaxation  in  this  frequency  range.  Increasing 
the  concent ration  increases  the  absorption  amplitude  and  results  in  either 
a  decreased  or  unchanged  relaxation  frequency.  This  lowering  of  the 

★ 

Tabulated  (or,  f)  values  are  available  from  G.  Atkinson  on  request. 
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relaxation  frequency  with  increasing  ion  concentration  makes  it  clear  that 
we  are  not  dealing  with  a  simple  bimoleeular  process.  The  individual  data 
points  represent  the  average  of  at  least  three  separate  determinations  of 
a  at  a  given  frequency  in  a  given  experimental  run,  and,  in  most  cases 
the  average  of  different  experimental  runs. 

Fig  2  shows  the  temperature  dependence  of  a  0.50  M  NagS04  solution. 

For  all  of  the  solutions  measured,  lowering  the  temperature  increased  the 
absorption  amplitude  and  decreased  the  relaxation  frequency.  The  shaded 
symbols  are  data  taken  with  the  high  frequency  cell  and  show  excellent 
agreement  in  their  overlap  with  the  medium  frequency  cell  data.  The 
relaxation  parameters  are  illustrated  for  the  25°  data.  The  data  are 
fitted  to  A,  f  ,  and  3'  using  a  non-linear  least  squares  analysis  programmed 
for  an  HP  9800  calculator.  In  no  case  aid  reach  a  constant  value 

at  high  frequencies.  However,  with  the  exception  of  the  5*5°  results,  (a'/fs 
did  reach  a  constant  low  frequency  value. 

Table  I  lists  the  relaxation  parameters  for  the  NagS04  solutions 
studied.  The  results  can  be  qualitatively  understood  if  we  assume  the 
existence  of  a  rapid  pre-equilibrium  for  the  formation  of  the  outer-sphere 
complex.  Then  the  relaxation  time  is  simply 


k5KIYl  ^CNa+  *  CS04~s^ 

1  +KIY+(CNa+  +  CS04'2) 


+  k_„ 


As  the  NagS04  concentration  is  increased  the  term  [k_  y+  (c^a+  +  cso  -2)] 


becomes  »  1  and 
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t'1  -  k2  +  k.?  (20) 

So  at  high  salt  concentration  t  becomes  independent  of  concentration. 

This  Is  seen  for  Na-jSO.1  at  25°  at  concentrations  greater  than  0.50  and 
at  15°  for  concentrations  greater  than  0. 50.  Although  it  is  useful  in 
explaining  the  above  observation,  eq  (19)  should  not  be  used  in  quanti¬ 
tatively  interpreting  our  data  since  it  is  based  on  the  assumption  that 
k.,,  k  »  k2,  k-r>.  This  is  unlikely  in  om*  case. 

The  changes  in  B'  with  concentration  and  temperature  can  be  related 
to  the  solvent  structuring  associated  with  "  structure- making"  and 
" structure- breaking"  ions^7  9  and  to  changes  in  the  excess  absorption 
due  to  the  higher  frequency  relaxation.  Na^SO-t,  classified  as  a  "structure 
maker, ,,2S  would  cause  increased  solvent  structuring  with  increased 
concentration  increasing  the  AB  value.  The  increased  association  at  the 
higher  concentrations  would  cause  the  A'  value  to  increase.  As  a  result 
we  would  expect  the  increased  B'  shown  in  Table  I.  A  t e t ramethy 1 armon i urn 
bromide  (TMA-Br)  gave  a  (cY'/fr)  value  that  was  independent  of  frequency 
indicating  no  measurable  association.  However,  a  1.5  M  solution  of  TMA-Br, 
the  same  ionic  strength  as  0.50  M  Na.^GO.jj  gave  a  B'  of  -5.5-  All  of  this 
is  attributable  to  the  solute- solvent  interaction  of  the  "structure 
breaking"  TMA-Br.*10  The  measurement  of  a  (0.60  M  TMA-Rr/O.  JOM  Ka,->S0.i) 
mixture  gave  values  of  (or'/f^')  that  are  3.0  to  1.2  units  lower  at  each 
frequency  than  the  pure  0.50  M  Na^SO-t  solution .  This  gives  A  and  f 
values  that  arc  the  same  as  the  pure  solution  but  a  B’  that  is  1.1  units 
lower.  For  the  same  mixture  at  15°.  (or'/f3)  values  are  lower  by  2.2  - 
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2.5  units  at  cadi  frequency.  tn  essence  the  added  TMA-Br  ban  a  drastic 
effect  on  AP  but  little  effect  on  A'  and  this  is  mirrored  by  the  changes 

in  B ' . 

Tn  order  to  compare  our  results  with  the  other  results  on  alkali 
metal  ions7'4  we  extended  the  ultrasonic  measurements  to  all  the  alkali 
metal  sulfates.1  Fig  3  shows  our  results  at  25°  and  0.A0  M  for  all  the 
salts.  As  in  Eigen's  work,  the  amplitude  of  the  sodium  salt  is  substantially 
greater  than  the  others  under  the  same  conditions.  The  relative  amplitude 
for  the  other  four  ions  depends  on  the  anion 

EDTA"4  Na  »  K  >  Li  ~  Rb  ~  Cs 

NTA“3  Na  »  Cs  ~  K  >  Li  >  Rb 

SO.!-2  Na  »  K  ~  Rb  ~  Cs  ~  Li  (This  work) 

Table  II  gives  the  fitted  relaxation  parameters  for  the  five  salts. 

The  results  show  one  important  difference  from  those  reported  by  Eigen  for 

the  large  multi  dentate  ligands.  Eigen  and  Mar.ss  found  relaxation  times  that 

increases  monatonically  from  Li1  to  Cs+  for  each  anion.  At  we  find 

this  same  order  of  increasing  relaxation  time  with  increasing  ionic  radius. 

however,  at  ?5°,  the  relaxation  time  for  the  I. i+  salt  has  gone  from  the 

slowest  to  the  fastest.  In  this  particular  ease,  the  three-parameter 

fit  yields  a  relaxation  frequency  of  5&0  MHz.  This  is  at  the  far  edge  of 

our  measurement  range  and,  therefore,  the  fitted  results  are  subject  to 

larger  errors  than  in  the  other  cases.  For  this  reason  we  also  tried  a 

two-parameter  fit  for  Li^SO*  fixing  B'.  For  a  B'  of  1.0,  A  has  decreased 

to  h.’j  and  f  to  505  MHz.  Wien  B'  is  increased  to  P.0,  A  becomes  3.6  and 

r 

f  decreases  to  397  MHz.  Since  the  measured  (&'/?')  is  *>  0  at  521  Mil  , 


However,  the  B'  of  Oji  obtained 


1U 

the  expected  D'  should  be  somewhat  Tower, 
with  the  three  -pa  rameter  fit  is  too  low.  The  best  value  for  the  relaxation 
frequency  of  Li^SO^  should  lie  between  the  r)60  'HI/,  value  and  the  339  MHz 
value  found  for  1.00  M  Li^SO^.  In  all  the  other  cases,  increasing  the 
concentration  either  lowered  the  f  or  left  it  unchanged. 

In  conclusion  we  feel  that  the  best  values  for  the  0.S0  M  Li^SO*  are 

those  from  the  two  parameter  fit:  B'  =  2.0  A  =  3.6  and  f  =  397  MHz.  This 

r 

still  makes  the  relaxation  time  for  LisSO^  substantially  faster  than  would 
be  expected  from  Eigen's  results.  Thus  at  25*  the  order  of  the  relaxation 
times  would  be 

Na  <  K  <  Rb  <  Cs  <  Li 

One  can  speculate  that  this  is  due  to  the  very  strong  hydration  of  the  Li 
ion.  This  diminishes  the  extent  of  inner-sphere  complexation  so  that  the 
bulk  of  the  absorption  results  from  the  much  faster  outer-sphere  process. 

The  multidentnte  ligands  used  in  the  Eigen  and  Maass  work  could  have 
reversed  this  effect.  At  5.5°  the  differences  caused  by  difference  in 
hydration  of  the  cations  is  diminished  since  the  entire  solvent  is  much 
more  structured.  This  permits  the  Li4  to  behave  in  a  more  "normal" 
fa shion . 

The  B'  values  for  the  K4 ,  Rb+ ,  and  Cs4  salts  arc  close  to  zero. 

This  is  consistent  with  our  qualitative  correlation  of  B'  with  the  structural 
effects  of  these  ions  and  the  contribution  of  the  fast  step  amplitude  to 
the  measured  amplitude.  All  three  of  these  cations  are  classed  as  "structure 
breakers"  while  S04  p  is  classed  ns  a  "structure-maker".  The  opposing 
cation-anion  effects  on  the  solvent  tend  to  cancel  each  other  leaving  only 
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a  small  amplitude  due  to  the  fast  step.  Li+  is  a  "structure-maker" 
so  we  would  expect  a  larger  B'  value  in  its  case. 

To  calculate  the  rate  constants,  we  will  use  the  equations  given  in 

the  theory  section.  First  an  'a'  value  must  be  chosen.  This  is  the 

"distance  of  closest  approach"  in  the  outer-sphere  complex.  This  is 

then  used  to  calculate  Ki.  and  k  .  Vie  have  assumed  k  is  indecendent  of 

I  “l  - 1 

ionic  strength  so  increases  in  Kj  are  reflected  by  increases  in  k| .  The 
Kj-  values  were  obtained  from  the  literature. 11  ’ 

Table  III  contains  the  most  important  results.  The  concentration 
dependence  of  the  rate  constants  for  HepSO*  reveals  the  same  unusual 
behavior  seen  in  the  relaxation  frequencies.  In  Table  II  we  saw  the  relaxation 
frequency  decrease  with  increasing  concentration.  Normally  no  dependence 
on  concentration  would  be  expected  since  this  is  a  uniraolecular  process. 

However,  in  this  system  the  coupling  between  k^  and  kp  is  closest  at  the 
lowest  concentration.  This  results  in  a  higher  kp  than  would  be  expected 
if  no  coupling  existed.  As  the  concentration  increases  k^  increases.  This 
tends  to  decouple  the  two  steps.  At  1.00  M  k^  and  kp  differ  by  a  factor 
of  16  so  the  two  steps  are  essentially  decoupled.  Calculation  of  k- 
for  1.00  K  Na2S04  using  the  rapid  pre-equilibrium  method  and  a  of  5-2  (ref  11) 
gives  5-1  x  10‘8  sec  1  instead  of  the  5-2  x  10a  sec  J  given  by  the  extended 
calculation.  Only  at  0.30  K  do  the  two  methods  give  appreciably  different 
results.  Vrhen  a  Kj.  of  6.6  is  used  (ref  35)  the  results  are  slightly  altered 
(see  Table  III). 

Since  the  value  for  ky  depends  on  the  value  of  'a'  assumed,  we 
varied  'a'  for  the  0.5  M  Ma  -SOn  calculations  at  PS°.  With  at  a  minimum 

O 

near  'a'  -  5  A  and  k  increasing  with  increasing  'a',  kp  is  smaller  at 
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both  4  and  6  A  than  it  is  at  5  A  .  If  one  '..•ere  to  suppose  that  'a'  did 
vary  with  concentration,  the  sine  would  be  inversely  related  to  the  con¬ 
centration.  The  effect  of  this  on  the  Ha^SO.j  results  would  be  to  give  a 
smaller  k2  at  lower  concentrations  while  leaving  the  high  concentration 
values  largely  unchanged.  This  would  bring  the  results  more  in  line  with  the 
more  aesthetically  pleasing  concentration  independent  rate  constant. 

To  examine  the  effect  of  temperature  on  k2  we  evaluated  the  results 

for  0.50  M  Na.oS04  at  15.0°  and  5-5°  using  'a'  =  5  A  .  The  values  for  Kj. 

at  these  temperatures  were  obtained  by  using  the  AH  value  of  +1.1  Kcal/mole 
reported  by  Austin  and  flair.32  The  calculated  k2  at  15.0°  is  very 
slightly  larger  than  the  value  at  25"  while  the  5*5°  value  is  appreciably 
less.  We  would  expect  the  decrease  of  kg  from  25"  to  5*5"  to  hold  true  for 
all  the  alkali  metal  sulfates  since  the  relaxation  frequencies  decrease 
for  all  and  all  should  have  positive  AH  values.  As  the  temperature  decreases 

k_A  also  decreases  due  primarily  to  the  decrease  in  the  diffusion  coefficients. 

The  two  parameter  fit  results  were  used  to  evaluate  the  rate  constants 
for  0.50  M  Li^SO-i.  The  Kj.  used  was  ^.5  (ref  55).  The  use  of  the  other  Re¬ 
value  in  the  literature  (5*9  in  ref  55)  altered  the  results  somewhat.  The 
results  are  given  in  Table  III.  The  initial  Kv  values  used  for  RbgS04  and 
Cs2S04  are  from  Reardon  (ref  55)-  However  it  was  necessary  to  use  a  K^. 
somewhat  larger  than  Reardon's  value  since  Reardon's  value  of  was  the 
same  as  our  K^. 

For  K2S04  two  different  Kj.  values  were  used.  The  value  of  9-1  (ref  ll) 
gives  k2  =  2.02  x  10°  sec  1  while  the  7-1  value  (ref  5*0  yields  k-  =  1.79  x 
109  sec-1.  Reardon's  work  indicates  that  the  smaller  value  of  is  a 
better  choice.  The  results  for  Rb2S04  er.d  CsgSO.j  are  very  similar  to  those 
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for  LioSO*  since  only  the  k  values  differ  appreciably.  This  effect  is 
due  to  the  larger  diffusion  coefficients  of  Hi/  and  Cs4 . 

Treatment  of  these  systems  using  the  simple  three-state  model  has 
greatly  facilitated  the  'understanding  and  analysis  of  the  data.  Fig  4 
gives  a  more  detailed  picture  of  the  Nn£S04  system  as  na  example.  However, 
pursuing  this  more  detailed  mechanism  in  our  calculations  does  not  seem 
fruitful.  In  essence  v:e  have  one  measured  relaxation  time  and  or.e  measured 
equilibrium  constant  for  each  salt  at  each  concentration  and  temperature.  In 
addition,  thermodynamic  data  is  lacking  for  the  secondary  equilibria  shown 
in  Fig  4.  In  all  of  these  systems  ve  are  dealing  with  electrolytes  that 
are  very  largely  dissociated.  In  addition  most  of  the  kinetic  processes 
are  close  to  diffusion  control.  Therefore,  we  do  not  believe  more  elaborate 
kinetic  analysis  is  fruitful  at.  this  time.  Giver,  the  experimental  concentra¬ 
tions,  the  assumptions  made  and  the  small  size  of  the  association,  constants, 
the  value  of  k2  at  0.50  M  for  Li,  Na,  Rb,  and  Cs  is  the  same,  1  x  10s  sec  1 . 
Potassium  is  higher  at  2  x  10H  sec  1 .  Since  all  but  the  K+  system  have 
similar  K^.  values,  this  is  not  a  surprising  result.  However,  it  is  surprising 
that  K2S04  with  a  larger  than  Na2S04  should  give  a  smaller  absorption 

amplitude.  The  only  explanation,  if  one  accepts  the  Ky  values,  is  that  the 

■f  "f* 

AV  for  the  K  association  is  much  smaller  than  for  the  Na  associat  ion. 

The  excess  absorption  values,  A,  can  be  used  to  calculate  the  values 
of  AV,  the  change  in  volume  on  ion  association.  Experimentally  we  obtain 
AVjj  the  volume  change  for  the  slower  relaxation.  AV,  is  not  available 
except  for  a  maximum  possible  value.  The  formation  of  the  outer-sphere 
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complex  should  have  n  .small  AV  since,  presumably,  no  water  is  displaced  in 

the  process.  Therefore,  the  contribution  of  this  coupled  step  to  £Mr,  the 

volume  change  for  the  second  step  will  be  small.  If  were  zero  or  t 

1  II 

<v~  Tjj  then  would  equal  AV^.  In  these  calculations  we  shall  always 
assume  that  the  contribution  of  AH  to  the  amplitudes  is  negligible. 

The  amplitude  expression  given  in  eq  7  can  then  be  simplified. 
Assuming  AH  -  0  and  using  the  thermodynamic  identities  36 


and 


CP  '  CV 


B 

K 


TV0*- 


Wliere  3  is  the  adiabatic  compressibility,  eq  7  becomes 

A  =  x  T 

RTVX 

Marks36  gives  values  for  the  velocity  end  density  of  the  alkali  sulfate 
solutions  studied.  X  can  be  evaluated  frees  the  known  concentrations  and 
equilibrium  constant  plus  the  calculated  activity  coefficients.  We  alrc 
assume  that  the  activity  coefficients  are  not  changing  appreciably  with 
concentration  in  this  range.  We  obtain  values  of  AVn  =  -4.4  cm ''/mole 
for  O.5O  M  NapSO^  at  25°  and  -6.6  cn3/rr.ole  at  5*5°.  The  for  1.00  K 

Na-504  at  25°  is  -5-3  cm3/r.ole.  Since  AV  changes  little  with  concentration, 
the  large  increases  in  A  with  increasing  concentration  must  bo  due  to 
increased  association  together  with  the  changes  in  solution  velocitv, 
density  and  relaxation  time.  Using  the  two  parameter  fit  for  0.50 
LioSO^  at  ?5°  we  obtain  AVjj.  =  -J.8  cm?/r.ole.  Tor  O.sO  M  K~S0.»  at  23® 
the  Kj.  of  9.1  gives  a  AVn  *  -4.4  cm3/nole.  The  Kv  of  7.1  gives  only  a 
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slightly  lower  value. 

To  obtain  the  best  AV2  we  should  do  a  normal  mode  analysis  of  AV-.  ^  to 
find  the  contribution  from  the  coupled  first  step?0  To  do  this  we  need  both 
eigenvalues  for  the  system.  These  can  be  obtained  by  using  the  rate  constants 
to  evaluate  the  coefficients  of  the  original  matrix  then  solving  for  the 
eigenvalues.  The  smaller  value  should  correspond  to  the  slower  step  and  the 
larger  one  to  the  rapid  first  step.  The  eigenvalues  are  then  used  to  find 
the  eigenvectors.  For  O.5O  M  IIa:JS04  at  25°  this  method  yields 


AV];  =  6.79  AVi  + 

1.00  av2 

AVII  =  0‘75A^  + 

1.00  AV2 

The  experimental  AV.^  is  -4.4  cm3/mole.  If  we  first  assume  AV^  =  0  then 
AVj  =  0.7  cm3/mole  and  AV2  =  -5*0  cm3/mole.  Even  if  AV^  is  allowed  to 
equal  all  the  residual  amplitude  (-7-2  cn3/mole)  then  AV*  =  -0.5  cm3/mole 
and  A V2  =  -4.2  cn3/nole.  So  even  a  large  change  in  A V_  has  a  relatively 
small  effect  on  our  AV2  value.  V.’e  could  confidently  expect  it  to  lie  between 
-4.2  and  -5-0  cm3/mole.  Thus  it  is  a  good  approximation  that  AV2  ~  AV^_, 
the  experimental  value,  and  that  the  two  steps  are  largely  decoupled. 

Furthermore  A V2  »  AVi  which  seems  eminently  sensible  when  we  think 
about  the  formation  of  inner-sphere  complexes  versus  outer- sphere  complexes. 
Table  IV  gives  the  AV  values  calculated  with  and  without  coupling  as  well  as 
showing  the  effect  of  different  AV^.  values  on  the  AVi  and  AV2  values.  It  is 
interesting  to  note  that  for  O.5O  M  Na2S04  at  5-5°>  £V2  =  -9-6  cm3/mole 
if  AVj  is  0  and  -8.5  cm^/nole  is  AVj  is  -7  cm'Vmole.  At  the  lower  temperature 
the  coupling  is  more  pronounced.  These  results  would  give  a  AV®  for  the  overal 
reaction  of  (AVi  +  AV2)  =  -  8.5  cm3/mole.  This  can  be  compared  with  the  much 
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larger  value  of  -15-8  cm3/mole  obtained  by  wester  and  Pytkowies3  for 
0.50  M  ha-SO*  at  2“  .  Their  value  was  obtained  from  the  pressure  dependence 
of  the  equilibrium  constant  determined  with  a  specific  ion  electrode.  One 
should  note  that  the  Kester  and  Pytkowica  AV°  for  I.ngSO*  is  even  larger  than 


that  determined  by  other  workers 


SB  >30  M-c 


for  I'.gSJ.j .  This  seems  improbable. 


Fisher  lias  measured  the  pressure  dependence  of  the  conductance  of  aqueous 
Ka2SO^  solutions.  His  recent  re-calculations41  gave  him  a  AV*  of  -5.25 
cra3/mole  for  Na+  -S04  2  ion  pair  formation.  This  is  in  very  good 
agreement  with  our  results.  Millero33  has  predicted  that  AV°  will 
increase  with  decreasing  temperature.  This  is  observed  in  our  work. 

Thus  a  substantial  amount  of  the  increased  sound  absorption  at  lower 
temperatures  is  due  to  an  increase  in  AV*  not  changes  in  density,  velocity 
and  degree  of  association.  These  last  three  effects  would  only  increase 
amplitude  by  about  15$  as  one  goes  from  25°  to  5.5" • 

In  conclusion  we  have  shown  that  alkali  metal  ion  association 
with  sulfate  proceeds  at  nearly  a  diffusion  controlled  rate.  The 
association  process  is  best  described  as  the  fast  formation  of  an  outer- 
sphere  complex  followed  by  rapid  conversion  to  inner-sphere  complexes. 

As  the  concentration  of  the  solution  decreases  the  rates  of  the  two 
steps  approach  each  other.  The  low  degree  of  association  gives  low 
amplitude  sound  absorption  with  most  of  the  amplitude  being  attributable 
to  the  outer-sphere  to  inner- sphere  conversion.  The  combination  of  low 
amplitudes  and  coupled  fast  steps  will  make  it  very  difficult  to  gain 
additional  information  about  these  systems  using  this  technique.  At 


higher  concentrations  it  does  appear  possible  to  treat  the  two  steps  as 
largely  decoupled  and  thereby  derive  very  useful  kinetic  information. 
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Table  I 

Relaxation  Parameters  for  Ha2S04  in  Water 


Concentration 

(Molarity) 

Temperature 

( °c) 

A 

fr(MHz) 

B' 

0.3CM 

25.0 

4.3 

313 

1.8 

0.5m 

25.0 

8.4 

221 

4.3 

0.70M 

25.0 

13.3 

221 

7-4 

1.0CM 

25.0 

22.7 

217 

9.8 

0.301 

15.0 

8.7 

293 

0.7 

O.5OM 

15.0 

15.9 

221 

4.3 

0.7CM 

15.0 

21.9 

153 

12.3 

1.00M 

15.0 

34.7 

155 

18.1 

0.30M 

5-5 

17.3 

233 

-0.7 

O.5OM 

5.5 

29.6 

169 

6.8 

Table  II 


Relaxation  Parameters  for  Alkali  Metal  Sulfates  in  Water 


Salt 

Concentration 

(Molarity) 

Temperature 

(*c) 

A 

fr(MHz) 

| 

B' 

Li2304 

0.50M 

25-0 

5.1 

560 

0.4 

LizSCU 

1.001 

25.0 

7-5 

550 

KrS04 

0.50M 

25.0 

5.0 

299 

1.6 

RbpS04 

O.5CM 

25.0 

7.2 

556 

-1.5 

Cs2S04 

0.50M 

25.0 

4.7 

520 

-0.5 

Li-S04 

0.50M 

5.5 

9.1 

142 

Li2S04 

1.00! 

5.5 

15.9 

128 

20.5 

K2S04 

0.50M 

5-5 

14.8 

178 

2.0 

Rb2S04 

0.50! 

5-5 

12.6 

161 

1.0 

C  s  -SO  4 

O.5OM 

5-5 

10.1 

219 

-1.1 

Rate  Constants  -Equilibrium  Constants  for  Selected  Alkali  Sulfates 
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Table  IV 

Reaction  Volumes  for  Selected  Alkali  Sulfates 


Salt 

Concen¬ 

tration 

(Molarity) 

Temper¬ 

ature 

C  c) 

AVII 

( cm3/mole) 

AVI 

(cra3/mole)  ( 

AVX 

cm^/mole)  ( 

&v2 

era '’/mole ) 

Na2S04 

0.5M 

25.0 

-U.h 

0 

+0.7 

-5.0 

-u.u 

-7.2 

-0.5 

-4.2 

Na2S04 

0.5M 

5-3 

-6.6 

0 

+1.7 

-9.6 

-6.6 

0 

t- 

1 

-0.1 

-8.5 

Na2S04 

1.0M 

25.O 

-5.3 

0 

+  0.5 

-10.6 

-5.3 

-10.0 

-1.5 

-9-6 

Li^S04 

0.5M 

25.O 

-3.8 

0 

+0.6 

-5.6 

CO 

K"\ 

1 

-5.2 

-0.2 

k2so4 

0.5M 

25.0 

-k.k 

0 

+0.5 

-2.5 

-4.4 

-6.0 

-0.2 

-2.1 

I 


Fig  1 


Fig  2 


Fig  3 

Fig  4 


The  excess  absorption,  a'/r?  x  101"  nepers  sec 2  era"1,  as  a 
function  of  the  log  of  the  frequency  for  sodium  sulfate  at 
25*C  and  varying  concentration. 

The  excess  absorption,  ff'/f2  x  lO1-*  nepers  sec2  cm  \  as  a 
function  of  the  log  of  the  frequency  for  0.50M  sodium  sulfate 
at  varying  temperatures.  The  values  from  a  three  parameter  fit 
at  25°C  are  A,  the  amplitude,  fr<  the  relaxation  frequency  and 
B',  the  residual  absorption  due  to  higher  frequency  relaxation 
and  soluto-^tolvcnt  interactions. 

Experimental  data  for  all  the  alkali  sulfates  at  25°C  and  0.50M 
ore  plotted  in  terms  of  o'/f2  x  lO1"’  nepers  sec2  cm  1 . 

Possible  additional  steps  for  the  ionic  association  of  sodium 


sulfate  are  shown.  The  indicates  an  outer  sphere  complex, 

while  the  inner  sphere  complex  is  represented  ns  a  contact  pair. 
The  charges  on  the  ions  have  been  omitted  for  simplicity. 


f  (MHz) 


THERMODYNAMICS  OF  IONIC  ASSOCIATION  IN  AQUEOUS  COPPER  SULFAT 
ION  SPECIFIC  ELECTRODE  MEASUREMENTS 

by 
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Al-Azhar  University 
Cairo,  Egypt 

and 

Gordon  Atkinson 
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Abstract 


The  stoichiometric  constants,  Kc,  for  CuS04  ion  association 

have  been  measured  at  15°,  25°,  and  55°  in  water,  aqueous  NaCl 

and  aqueous  NaC104  using  a  divalent  ion  electrode.  The  K  values 

c 

were  converted  to  the  infinite  dilution  values  using  standard 
activity  coefficient  methods.  The  values  found  by  this  method 
are  significantly  higher  than  those  measured  spectrophotometrically. 
This  difference  is  attributed  to  the  fact  that  the  electrode 
technique  can  be  used  to  measure  both  inner  and  outer  sphere  ion 
pairs  while  the  spectral  method  only  measures  the  inner  sphere  pairs. 
Combining  the  results  from  the  two  methods  allows  us  to  calculate 
the  thermodynamic  parameters  for  the  individual  steps.  At  25° 
we  find  =  99-2  (free  ions  to  outer  sphere)  and  =  1.62 
(outer  to  inner  sphere) . 


Int  ro  l  uct  lor. 


Recent  years  Have  seen  a  greatly  increased  interest  in  ion 

specific  electrodes.  The  utility  of  the  now  classical  glass 

electrode  selective  for  hydros. ium  ions  was  demonstrated  in  the 

late  1950’s.1  The  development  of  glass  electrodes  specific  for 

^  * 

univalent  cations  occurred  in  the  late  19?0’s.  ’  But  the 

development  in  the  19o0’s  of  liquid  membrane  and  solid  state 

electrodes  for  a  wide  variety  of  cations  and  anions  marked  the 

k 

beginning  of  the  new  era.  The  new  electrodes  together  with 
high  precision  solid  state,  digital  pH-mv  meters  makes  broader 
application  of  potentiometry  in  physical  and  inorganic  chemistry 
a  certainty. 

We  have  been  concerned  with  the  measurement  of  ion  associa¬ 
tion  for  years.  Although  numerous  techniques  have  been  used  to 
measure  association  in  classical  systems  such  as  MgS04  and  CuS0.» , 
the  situation  is  still  not  very  satisfactory.  Cne  of  the  problems 
is  the  tendency  of  investigators  to  choose  a  single  experimental 
technique.  This  makes  it  very  difficult  at  times  to  compare 
results  on  the  same  system  done  by  different  investigators.  The 
differences  in  ion  association  constants  for  CuSO.»  in  water  at  25* 
seem  large  compared  to  the  precision  of  the  different  methods. 

Therefore,  we  decided  to  apply  a  variety  of  methods  to  the  same 

/ 

system.  In  the  first  paper  "  we  repeated  and  extended  the  use 
of  UV-visible  spectroscopy  to  measure  ion  association  in  aqueous 
CUSO4.  In  this  paper  we  will  examine  the  same  system  using  the 
poter.tiometric  technique  with  a  ion  specific  electrode. 


Experimental 


A .  Materials 

1.  C113O4:  This  was  Fisher  Reagent  Grade  recrystallised 
twice  from  distilled  and  deionised  water. 

2.  CuBDS  (Copper  (ll)  m-tencer.edi  sulfonate^ 

Crude  H2BDS,  as  obtained  from  Eastman  Kodak,  contains 
excess  H2SO4.  The  crude  H2BDS  in  aqueous  solution  is 
treated  with  Ba(OH)2  until  no  more  5aS04  is  formed.  The 
resulting  BaBDS  solution  is  passed  through  a  cation  exchange 
resin  in  the  hydrogen  form  to  obtain  a  pure  H2BDS  solution 
which  can  be  standardised. 

To  prepare  CuBDS  equivalent  amounts  of  the  pure  acid 
solution  and  dried  CuC03  are  reacted.  The  resulting 
solution  is  evaporated  until  crystals  form.  The  product 
is  recrystallised  twice  from  deionised  water.  Purity  is 
checked  by  EPTA  analysis  for  the  Cu  and  ion  exchange  for 
the  BDS . 

5.  NaCl:  This  was  Fisher  Reagent  grade  dried  at  105°C. 

NadC^:  Fisher  Reagent  grade. 

B.  Measurements 

The  potentiometric  measurements  were  made  using  an 
Orion  Model  801  pH-nv  meter  equipped  with  an  Orion  Model 
92-32  divalent  cation  electrode  and  a  90-01  single  junction 
reference  electrode.  The  92-32  was  modified  for  Cu*1- 
selectivity  using  the  internal  and  filling  solutions  suggested 
and  supplied  by  Orion. 

Measurements  were  carried  out  in  a  large  volume  jacketed 
glass  cell  whose  working  solution  temperature  could  be 
controlled  to  +0.05°. 


All  working  ar.d  reference  solutions  were  preoared  by 

4-0 

weight  dilution  techniques.  The  Cu  **  concentration  of 
each  solution  was  checked  by  EX'A  titration. 


L. 


The  primary  problem  in  ion  specific  electrode  measurements 
is  establishing  a  standard  curve  of  potential  versus  activity.  We 
elected  to  use  Cu3DS  as  our  standard  reference  because  of  the 

rj  Q 

strong  evidence  ’  that  it  is  essentially  completely  unassociated 
in  aqueous  solution.  It  is  therefore  advantageous  to  use  it  as  an 
activity  reference  for  CUSO4  since  it  is  also  a  2-2  salt.  The  use 
of  2-1  salts  such  as  Cu(N03)2  as  reference  presents  additional 
problems  in  activity-concentration  conversions. 

A  standard  reference  curve  was  run  before  and  after  each 
set  of  ''unknown''  measurements.  Unless  the  curves  agreed  the 
data  were  discarded.  It  is  our  experience  that  the  liquid  mem¬ 
brane  electrode  cannot  be  trusted  to  be  stable  for  more  than  1-2 
hours. 

Three  different  media  were  used  in  the  CuS04  measurements: 
pure  water,  aqueous  NaCl,  and  aqueous  NaC104.  Standard  reference 
curves  were  run  with  Cu3DS  in  the  same  three  media.  All  systems 
were  examined  at  15°,  25°,  55°. 

The  activity  of  Cu+2  in  the  ''unknown''  solutions  was  cal¬ 
culated  from  the  standard  reference  curve  and  converted  to  [Cu+2] 

9 

using  the  known  activity  coefficients  of  CuS04  at  the  given  ionic 
strength.  The  (  [cu+2],  I  )  values  must,  of  course,  be  obtained 
by  an  iterative  procedure. 

All  measurements  were  made  at  a  pH  of  4  to  remove  the  necessity 

of  correcting  for  Cu+2  hydrolysis. 

Table  I  is  an  example  of  a  data  set  showing  the  calculation 

of  the  stoichiometric  association  constants,  K£,  for  CuS04  in  H20 

at  25°C.  Table  II  summarizes  the  K  results  in  the  three  media. 

c 

Evaluation  of  Association  Constants 

It  is  now  useful  to  convert  the  stoichiometric  constants,  Kc, 
to  the  infinite  dilution  state  values,  K^. 


Table  II 


CuS04 


Cu-  SOi 


[CuS04]  vp 
[Cu2*]  [S042']  y; 


K  /y2 
c'  + 


assuming  y  =  1 


We  have  used  the  some  methods  described  in  paper  I.  In  the  first 

method  the  Davies  equation10  is  used  to  calculate  y  Dor  each  K 

-  c  ] 

at  a  given  I  and  an  average  is  ootained.  In  the  second  method 
eq  (l)  is  rearranged  to 


log  Kc  -  log  y;  =  log  Ka 


Then  the  extended  Debye-Huckel  equation  L  is  used  to  calculate  y+ 
at  different  I  values  with  a  series  of  a0 ,  mean  distance  of 
closest  approach  values.  Then  a  plot  of  the  L.H.S.  of  eq(2)  versus 
I  gives  log  as  the  intercept.  Both  techniques  can  be  made  to 
agree  to  within  5^-  The  values  and  the  other  thermodynamic 
parameters  are  given  in  Table  III  (pure  H20). 

It  is  interesting  to  compare  the  values  in  the  three  different 


media  at  25* 


Medium 


Na?C104  (aq) 
NaCl  (aq) 


ka  (c_1) 
259-8  +  11.5 

252.5  +  12.8 


21U.U  +  18.2 


Although  the  values  are  close  to  being  the  same  within  experimental 
error,  they  are  somewhat  different.  These  apparent  differences 
could  arise  in  two  general  ways.  First,  the  K^'s  have  been  cal¬ 
culated  using  a  pure  electrolyte  activity  coefficient  approach, 
that  is  assuming  y+  depends  only  on  ionic  strength.  This  is  known 
to  be  generally  not  true.  However,  appropriate  corrections  would 


be  difficult  to  make  considering  the  overall  precision  of  this 

data.  The  second  possible  cause  for  the  differences  is  our  neglect 

13 

of  other  weak  ion  association  equilibria  in  particular, 


ar.d 


Na1+(aq)  +  S042"  (aq) 
Cu2’*  (aq)  +  Cl"  (aq) 


_ \ 

_ \ 

^ - 


t>aSO|]  (aq) 
[cuCl+]  (aq) 


In  the  NaC104  media,  the  first  equilibrium  would  decrease  the 
S04~  available  for  Cu2^  association.  In  the  NaCl  media,  both 
equilibria  would  effectively  remove  reaction  partners  in  the 
Cuph-  S04  association  process.  Semi-quantitative  corrections 
made  at  various  Kc  values  in  the  two  media  are  of  the  right  order. 


Discussion 

It  is  now  appropriate  to  compare  the  values  for  [cuS04] 
ion  pair  formation  obtained  by  the  I.S.E.  technique  with  these 
previously  obtained  spectrophotometrically.  At  25°  we  find 


C  =  160.6  C  1  ( spectrophotometric) 

\  Ka  =  259.8  C"1  (I.S.E.) 

lU 

In  paper  I  it  was  suggested  that  this  could  be  explained  qual¬ 
itatively  by  the  fact  that  [cu^  -  S042"]  association  takes  place 
in  a  step-wise  fashion1'.  In  the  simplest  form  this  means 


Step  j 

Cu2*  (aq)  +  SOf”  (aq)  ^  [cu^ - SOf"] 

STATE  ©  © 


II 


[cursor] 

0 


where  state  © 
state  @ 
state  0 


free  ions 

outer-sphere  ion  pair 
inner-sphere  ion  pair 


The  comparison  of  the  various  association  concepts  is  conceptually 
simplified  by  considering  a  pure  CuS04  solution.  If  we  let 


a.^  =  Activity  of  species  i 

then  K  = 


a- 


ai 


II 


ap 


Now  the  "real"  thermodynamic  must  be  for  the  formation  of  an 
equilibrium  mixture  of  ©  and  from  0 


K. 


a2  +  a3 
a? 


KI  +  KI  KII 

This  should  be  equivalent  to  the  KA  determined  by  I.S.E.  since  in 
this  method  we  measure  ai  and  obtain  the  sum  of  (  a2  +  a3  )  from 
stoichiometry. 

However,  in  the  spectrophotometric  method  we  measure  C3 
(  corresponding  to  a3  )  directly  and  assume  everything  else  is 
"free"  ion.  Therefore,  this  method  lumps  ^  and  ^together. 


So  Ks  *  ka 


( spectrophotometric) 


a3 

(ai  +  a 2) 2 


We  can  now  show  that 


Ks  -  KIKII 


The 


Therefore, knowing  both  K.  and  K_,  we  car.  calculate  KT  and  KTT. 

Pi  O  1  II 

results  are  given  in  Table  IV . 


Table  IV 

T  (°C) 

Ks  (C1) 

ka  (C_1) 

Kj  (C'1) 

KII 

15 

145.0 

235. 

90.4 

1.604 

25 

160.6 

259.8 

99-2 

1.619 

35 

169.6 

274.2 

104.6 

1.621 

45 

201.4 

_ 

“ 

- 

With  real  trepidation  because  of  the  large  errors  involved,  we 
can  use  these  values  to  estimate  the  AH  parameters  and  other  para¬ 
meters  of  the  two  steps. 

Note  that 


AX» 

=  AX°  +  AX° 

I  II 

while 

AX; 

KII 

=  AX°  +  11 

AX° 

II 

for  any  parameter  X 


The  derived  25°  values  are  given  in  Table  V. 


The  AX  parameters  are  reasonably  consistent  with  those  calculable 
from  the  K.  and  Kg  data,  but  clearly  cannot  be  trusted  very  far. 
Neither  experimental  technique  is  very  high  precision.  However, 
the  results  do  encourage  our  belief  that  the  different  experimental 
techniques  do,  indeed,  measure  different  aspects  of  ion  association. 

Another  interesting  comparison  can  be  made  if  we  compare  the 
results  from  this  work  with  results  from  previous  ultrasonic 
investigations  of  step-wise  association  in  2-2  sulfates  .  If 
we  use  MgS04  as  a  contrast,  we  can  lump  the  first  two  steps  of 
this  three-step  system  together  to  obtain  a  valid  comparison  with 
CuS04  where  only  two  steps  have  ever  teen  observed. 

CuS04  MgS04 

Free  Ion  to  Outer  Sphere  Conversion  K  99.2  15^.2 

Outer  to  Inner  Sphere  Conversion  K  1.62  0.17s 

So  we  find  a  much  greater  conversion  of  outer  sphere  to  inner  sphere 
ion  pairs  in  the  case  of  CuS04.  This  can  be  related  to  the  very 
high  rate  of  H20  exchange  on  Cu2*  as  compared  to  Mg21". 

In  conclusion,  we  have  shown  that  ion  specific  electrodes 
can  be  used  for  reasonable  precise  determinations  for  the  thermo¬ 
dynamics  of  ion  association  of  divalent  metals.  The  method  has 
an  inherent  advantage  in  that  it  determines  free  ion  activity 
ignoring  all  association  processes.  A  comparison  of  I.S.E.  values 
with  these  obtained  by  a  method  sensitive  only  to  inner-sphere 
complexation  can  be  used  to  obtain  useful  step-wise  association 
parameters. 
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Introduction 

The  problem  of  ion  association  in  electrolytes  has  been  an  important 
subject  of  investigation  for  physical  chemists  for  more  than  l60  years.1  3 

Theories  of  ionic  association  began  when  Volta  discovered  a  way  of 
producing  electricity  by  means  of  chemical  reactions.  After  that,  Faraday, 
Clausius,  Arrhenius  and  Ostvald  made  substantial  experimental  studies  in  the 
last  century.  However,  the  first  mathematical  approaches  to  the  problem 
were  begun  in  this  century  by  Bjerrum,4  Muller,5  Gronwell,6  Fuoss"'  and  Denison 
and  Ramsey.8  Detailed  explanations  of  the  full  mathematical  derivations  of 
the  various  models  can  be  found  elsewhere.9 

If  we  turn  to  the  experimental  methods  by  which  one  can  study  ion 
association,  several  ^  have  been  in  common  use.  The  most  important  have 
been  conductance,  potenticmetric,  spectrophotometric,  polarographic,  ion- 
exchange,  solubility,  solvent  extraction,  Raman  spectra  and  reaction  kinetics. 

It  is  desirable  to  study  the  same  system^  by  as  many  techniques  as  possible^ 
so  that  the  investigator  can  make  meaningful  comparisons. 

Copper  sulphate  has  been  studied  by  many  authors^  using  several  techniques 
including  conductance,  E.M.F.  and  spectroscopic.  However,  most  of  the  different 
measurements  have  been  carried  out  under  different  conditions  so  that  comparison 
becomes  difficult  and  ambiguous.  We  shall  only  review  the  results  obtained 
using  the  two  major  methods,  conductance  and  spectrophotometric. 

The  thermodynamic  association  constant  (Ka)  of  CuS04  in  water  at  ?5''C 

obtained  from  conductance10'11  measurements  ranged  from  200  to  232  M  1 .  The 

G*/W5.  o-i 

spectrophotometric  method  rei” i'  ■.  values  for  ranging  from  125  to  2o5  M  . 


1 


obtained  values  of  212.7 


and  222.2 


For  example,  Monk  and  co-workers  *■ 
vhile  Nasanen13  obtained  a  value  of  125  M  1 .  The  latter  author  could  not  find 
an  explanation  for  the  large  difference  between  the  results. 

Prue1-'  obtained  a  value  of  285. 7\.  while  Matheson15  could  not  obtain  a 

V 

V 

unique  value  for  without  arbitrary  conditions.  He  obtained  values  in  the  re 
of  170-250  M  1  •  Later,  Petrucci  and  KeTm.es  17  made  LTV  spectroscopic  studies 


on  CuS04  at  25°C  and  obtained  a  value  of  220-2J0  M  1  for  K, . 

t i 


It  is  obvious  that  the  values  obtained,  even  using  the  same  experimental 


method,  differ  greatly  from  author  to  author.  It  is  also  noticeable  that  the 
values  of  K,  obtained  by  the  soectroohotometric  method  are  smaller  than  those 
obtained  from  the  conductance  method. 


This  paper  and  the  one  which  follows  are  the  result  of  thorough  investi¬ 
gations  in  our  laboratories  on  CuS04  in  various  aqueous  media  and  under  a  wide 
range  of  experimental  conditions.  In  this  communication  the  spectrophotometric 
measurements  are  reported  at  four  temperatures  and  for  three  ionic  strengths 
at  each  temperature.  The  following  p>aper  includes  the  ior.-selective  electrode 


measurements  and  a  comparison  of  the  results  found  using  both  techniques. 


Experimental 


A .  Materials: 

Copper  perchlorate  was  prepared  from  ar.alar  ccpper(II^oxide  using  a  sli 
excess  of  analar  perchloric  acid,  then  filtering  and  crystallising.  The  produ 
was  then  recrystallised  from  0.01  M  KCIO4 .  Sodium  sulphate  was  Baker  analytic 
reagent.  Sodium  chloride  was  Fisher  reagent. 

B .  General  Equipment : 

Measurements  were  made  using  a  Zeiss  QM3  spectrophotometer  equipped  wit 
a  1  cm  cell.  The  temperature  of  the  solution  was  maintained  constant  by 
circulating  water  through  the  cell  compartment.  The  temperature  was  cor.trolle 

to  +  0.01°C. 


C.  Measurements  and  Collect ir.g  Bata: 


Optical  density  (D)  was  determined  at  each  temperature  using  1  cm  cell 
for  the  solution  ( D ’ )  and  the  reference  (D).  Copper  perchlorate  was  used  as 
the  reference  and  the  sample  contained  Na-'S04 ,  Cu(C104)-  and  N aC  10 4 •  (D-D') 


was  read  directly  and  the  measurements  were  repeated  several  times. 


Bupli cat  e 


solutions  of  s ancle  were  measured  and  the  values  reported  are  the  average  of  f 


two  readings  for  each  concentration. 

To  check  any  possible  error  in  the  measurements ,  a  series  of  Cu(ClC.;\- 


solutions  were  carried  out  at  230  nn  wavelength  using  an  assortment  of  cells  o 
different  path  lengths.  The  results  showed  that  the  optical  density  over  the 
range  of  this  work  is  good  to  within  l*c.  The  source  used  for  the  UY  work  was 


3 


1. 


the  hydrogen  lamp. 

The  cell  used  was  1  cm  length  quart;:  and  was  soaked  in  HSO*  over 
night,  cleaned,  ami  washed  with  distilled  water,  acetone,  then  dried  with  air 
before  use.  Stock  solution  of  Cu(ClO.i).-  was  prepared  from  reagent  grade.  It 
contained  the  required  amount  of  HC10*  to  prevent  the  hydrolysis  of  Cu. 

Seme  exploratory  experiments  were  carried  out  to  study  the  possible 
hydrolysis  of  Cu  •  ion  and  determine  the  absorption  peak.  The  results  are 
given  in  Table  I  as  the  absorbance  of  Cu(C10.»).->  in  presence  of  varying  amounts  of 
HCIO4  at  25#C.  It  is  seen  that,  the  absorbance  does  not  change  after  the 
concentration  of  HCIO4  reaches  a  value  above  1 . ft  x  10  '  M.  This  means  that  all 
our  measurements  must  he  carried  out  in  the  presence  of  at  least  this  amount  of 
HC10.t  to  prevent  hydrolysis  of  Cu^‘  at  05"C.  In  all  the  experiments  HCIO4 
was  added  to  minimise  this  hydrolysis. 


Results  and  Discussion 


The  primary  aim  of  this  work  is  to  study  the  effect  of  temperature  and 
medium  on  the  association  constant  determined  spectrophotometrically.  CuSO* 
has  a  particular  interest  both  because  it  is  a  classical  case  of  D-2  electrolyte 
association  and  because  the  association  can  be  measured  directly  by  the  spectro- 
photometric  technique.  Secondly,  it  is  necessary  to  obtain  more  precise  data 
under  varying  conditions  so  that  the  apparent  differences  in  determined  by 
the  different  methods  can  be  examined  in  detail.  Finally,  the  present  work  is 
intended  to  lay  a  foundation  for  the  investigation  of  the  same  system  at 
high  pressures. 

As  mentioned  earlier,  the  method  consists  of  measuring  the  difference  in 
optical  density  between  two  solutions,  a  reference  containing  Cu(C10^).'  - 
NaC104  and  a  jubrnMe  solution  of  Cu(C104)r  -  N a rS0 4  and  enough  NaClO«  to  mainiain 

'V  Lsl'T-  <>»• 

constant  ionic  strength.  a  is  the  [Cu(C104) sto i oh J  —  sloi  ch ' 

•  V  (.fV 


c  [l!aC104] 


r.toich 


rCuS04]  (ion  pair).  The  reaction  under  investigatior 


is  represented  hy  the  equatft 


Cu*2  +  S04‘ 


Concentration  (a  -  x)  (b  -  x) 


sCuS04 


It  is  well  known  that  both  free  Cu  2  and  thtKCuS04  ion-pair  absorb  in  the  I'V 
range.  The  measurements  were  carried  out  at  a  wavelength  of  r'0  nm.  If  the 
Cu*'  (aq)  and  rCuFCht]  arc  the  only  two  species  absorlnSg  at  Dc'0  rm.,  then  using 
Beer's  law  we  have  D  eln  but  1  --  1  ;  then  D  ea  ;  wheiS^D  is  the  optical 
density  of  Cu*  ',  e.  is  the  molar  extinction  coefficient  of  Cu  '  (iiqV>«^  is  the 
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concent  rut  ioH^of  Cu  ‘’(aq)  and  D'  =  e'x  +  e  (a  -  x)  vj>^re  I)'  =  optical  density 
of  [CuS04]  and  tnK^Cu’ ;’( aq)  ions,  e'  =  extirpation  coefficient  of  [CufO.,], 
a  -  x  =  [Cu+~]  in  the\anple  solution^-D'  =  Ae  •  x.  But  the  stoichiometric 
association  constant  of  CulB^i^given  by  =  x/(a  -  x)  (b  -  x). 

Combining  the  above  equjvj^fon,  w^obtain 


2c  =  a  +b  -  x 
-  D’  Ae 


Ae  • 


where  Ae  -—-'e  -  e  ' 

Our  task  now  is  to  calculate  both  and  Ae ,  since  the  other  parameters 
in  the  last  equation  are  known  (except  x).  Two  methods  were  used  to  calculate 
K  and  Ae .  The  first  is  the  graphical  method,  in  which  x  is  neglected  first 
and  ab/D-D'  is  plotted  against  (a  +  b) .  From  the  slope  and  intercept  both 
Ae  and  Kc  are  obtained.  A  value  of  x  is  the  obtained  by 

x  =  (a  *  b)  +  -  -D— D,  *  Ae 

c 

then  a  new  graph  is  plotted  between  ab/D-D'  against  (a  +  b  -  x).  New  values 
of  K^,  and  Ae  are  obtained.  This  process  is  repeated  several  times  until 
constant  values  of  x,  and  AC  are  obtained.  Three  iterations  were  necessary 
to  obtain  such  constant  values. 

The  second  method  in  calculating  and  Ae  is  simply  using  a  least 
squares  analysis  with  an  iterative  procedure  similar  to  that  used  in  the 
graphical  method.  A  computer  program  was  written  and  used  in  this  case.  The 
results  obtained  from  both  methods  are  essentially  identical.  Tables  II  and  III 
give  examples  of  the  data  analysis  at  ?2°C  and  respectively,  while  Fig. 

1  shows  the  plots  of  ab/D-D'  against  (a  +  b  -  x)  at  ?5°C  for  three  .ionic 
strengths.  Table  TV  reoresents  the  stoichiometric  association  constants 


(K)  at  various  ionic  strengths  and  temperatures.  To  correct  for  the  ionic 

V* 
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strength  effects  the  stoichiometric  association  constants  at  the  various 

ionic  strengths  must  be  extrapolated  to  zero  ionic  strength.  This  simply 

means  calculation  of  the  thermodynamic  association  constants  (K.)  from  K_ 

A  C 

values.  To  carry  out  such  a  calculation  activity  coefficient  values  must 


be  known  or  estimated.  is  related  to  Kj,  by  the  following  relationship 


KA  =  'V 


which  can  be  rewritten  as 


CuSCU 


Xu2  • 


ka  =  KC  '  -f" 

A  c  r-2  V*. 

where  ^  _  =  activity  coefficient  of  [CuS04]  ion  pair  is  considered  to  be  unity 

Yc,+2  *»  ^0  =  Xl  the  mean  ionic  activity  coefficient  of  the  free  ions. 

In  our  calculations  for  obtaining  from  the  mean  ionic  activity  coefficients 
for  the  ionic  species  were  obtained  using  the  Davies  equation 

-  Log  %  =  0.5  Z?  J  ,  -  0.3l! 

C  1  +  r 

where  Z^  =  charge  of  ion  and  I  is  the  ionic  strength.  The  temperature  effects 
on  the  activity  coefficients  are  not  considered.  The  values  in  Harnea  and 
Owen's  tables  indicate  that  for  similar  salts  varies  less  than  1.5$  in  a 
30oC  range.  The  mean  activity  coefficients  of  the  Cu+,?(aq)  ions  as  a  function 
of  temperature  are  notNavailable. 

is  calculated  at  the  same  temperature,  using  all  the  ionic  strengths 
at  this  temperature  and  an  average  value  is  obtained  along  with  the  standard 


error. 

t 

Another  way  of  obtaining  the  thermodynamic  association  constant  was 
developed  by  Petrucci.  Here  a  plot  of  (Log  K  -  log  Y_S)  versus  the  ionic 
strength  is  made  according  to  the  equation 
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Log  Kc  -  log  &2  =  Log  Ka 

In  this  method  several  different  values  of  are  calculated  by  the  extended 
Debye-Huckel  activity  coefficient  equation  using  different  values  for  closest 
distance  of  approaches0.  When  three  values  for  a°  were  used(o,  2.5  and  5  A^ 

A  6<f- 

they  gave  a  thermodynamic  association  constant  within  5^  f-ee«  the  value 

obtained  by  the  other  method.  Table  V  shows  the  values  at  15,  25,  32  and 

l*5°C.  Also  included  in  the  same  table  are  the  thermodynamic  parameters  £H° , 

/vi °  and  AS°  at  the  corresponding  temperatures. 

Then  comparing  these  results  with  other  workers  using  the  same  method  we 

find  that  the  agreement  is  very  good  with  Petrucci's  results  at  25°C  for  the 

same  ionic  strengths  (0.06  and  0.1).  Table  VI  summarizes  this  comparison. 

However,  the  values  obtained  spectrophotometrically  are  distinctly 

lower  than  the  values  obtained  by  the  other  methods.  Qualitatively  this  can  be 

rationalized  by  considering  that  CuS04  ion  association  is  a  multi-step 

process  whose  most  simple  adequate  representation  is 

KI  KII 

Cu+2(aq)  +  S04"2(aq)  ^  [Cu+2.  .  . S04'2]  ^  [Cu2S04'2] 

State  ,  0  , 


Here  State  1  =  free  ions 

State  2  =  outer-sphere  ion  pair 
State  5  =  inner-sphere  ion  pair 

Such  behavior  has  been  amply  demonstrated17  for  many  such  systems  using  the 
methods  of  relaxation  kinetics.  Therefore,  in  such  systems  there  is  an 


I 
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equilibrium  mixture  of  free  ions,  outer-sphere  ion  pairs  and  inner-sphere 
ion  pairs.  Then  in  comparing  the  results  of  different  experimental  methods 
for  it  is  important  that  we  recognize  what  species  the  technique  measures. 
Basically^methods  such  as  conductance  and  E.M.F.  measure  the  free  ion^  concen¬ 
tration  or  activity.  The  ion  pair  concentrations  are  then  obtained  by  difference 
(and  by  theory!).  However,  the  spectrophotometric  method^particularly  in  the 
UV^can  distinguish  between  inner-  and  outer-sphere  ion  pairs.  Smithson  and 
Williams10  pointed  out  in  1958  that  outer-sphere  ion  pairs  would  probably  have 
little  effect  on  the  UV  spectrum  of  the  absorbing  species. 

In  the  case  of  [CuS04],  then,  the  spectrophotometric  method  will 
"see"  the  outer-sphere  ion  pairs  as  free  ions.  This  will  lead  to  a  lower 
value.  In  principle  the  combination  of  this  KA  value  and  a  KA  from  a 
conductance  or  E.M.F.  method  would  allow  the  calculation  of  and  step¬ 
wise  association  constants  and  the  relative  populations  of  inner-  and  outer- 
sphere  species.  This  technique  will  be  elaborated  on  in  the  following  paper. 
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Table  II 


Stoichiometric  Association  Constant  and  Spect rophotomet ri c  Data 
of  CuSC4  at  32°  C  [Cu(C104 ) -»]  =  a  =  4.05  x  10~3  M 
[HCIO4]  =  1.01  x  10" 3  M 


[NaaSO*] 

M  x  103 

(b) 

[NaClO-t] 

M  x  103 

(c) 

D-D' 

ab 

(a  +  b  -  x) 

M  x  10® 

D-D' 

x  10  4 

4 

67.8 

0.092 

1.76 

8.1 

3 

56.5 

0.167 

1.94 

12.97 

12 

45.2 

0.237 

2.0S 

15.9s* 

16 

>3.9 

0.298 

2.17 

19.19 

20 

22.6 

0.237 

2.33 

23.51 

24 

11.3 

0.594 

2.46 

27.02 

23 

0.0 

0.438 

2.59 

50.54 

^  =  0.590; 

Ae  =  270.3; 

k  =  25.34  M'1; 

I  =  0.091  +  0.001 

Table  III 


Stoichiometric  Association  Constants  and  Spectrophotometric  Data 
of  CuS04  at  1*5 °C  [Cu(C104)r]  =  a  =  5.0U  x  10' 3M 
[HC104J  =  7.587  x  10"“  M 


[Na2S04] 

M  x  103 

(b) 

[NaC104] 

M  x  10 3 

(c) 

D-D' 

ab 

(a  +  b  -  x) 

M  x  103 

b  ‘ii 

> — t 

1 

X 

Q 

5 

25.6 

0.119 

7-66 

6.1*12 

U 

21.0 

0.165 

7-46 

5.722 

6 

15-7 

0.216 

8.44 

9-101 

8 

10.5 

0.277 

8.79 

10.508 

10 

5-2 

0.518 

9.56 

12.965 

12 

0.0 

0.564 

10.02 

ll*.550 

y  =  0.185; 

0) 

CO 

bT\ 

II 

W 

<3 

K  =  50.8b  m'1} 

I  =  0.0U18  + 

0.002 

Table  IV 


Stoichiometric  Association  Constants  of  CuSO^ 


at  Various  Ionic  Strengths 


and  Temperatures  from  Speccrophotcmetric  Methods 


Temp. 

Ionic  Strength 

Ac 

KC 

°C 

I 

lm"1 

15 

0.0435  +  C.OOUl 

215.1 

51-63 

0.0634  t  °-°°15 

200.0 

27.59 

0.0Q22  +  0.0025 

200.0 

22.95 

25 

0.0421  +  0.002T 

250.0 

55-10 

0.0655  t  °-005° 

209.5 

51.97 

0.0910  +  O.OO35 

242.1 

24.43 

52 

0.0429  +  0.0024 

266.60 

38.26 

0.0627  +  0.0023 

234.15 

55-55 

0.0912  +  0.0010 

270.27 

25.54 

45 

0.0418  +  0.0020 

558.93 

50.86 

0.0620  +  0.0019 

248.43 

56.79 

0.0995  1  0.0017 

384.61 

28.57 

The  errors 

in  Ae  and  are  +  5'«. 

Additional 

digits  are  carried  in  the 

results  to 

eliminate  rounding 

off  errors  in  other  calculations  . 


Table  V 


Thermodynamic  Parameters 

of  Ion  Pair 

Formation 

of  CuS04 

T 

ka 

£H° 

£G° 

AS° 

°C 

Uf1 

Kcal  M'1 

Kcal  M~ 1 

cal  M~ 1  K° 

15 

145.0 

2.U0 

-2.86 

18.0 

25 

160.6 

2.  hO 

-5.02 

18.0 

32 

169.6 

2.U0 

-3-12 

18.0 

^5 


201.  b 


2. hO 


-3-39 


18.0 


Table  VI 


Summary  of  the  Spectrophotometrically  Determined  Stoichiometric 
Association  Constants  of  CuS04  at  25°C 
in  Aqueous  NaC104  Media 


I 

Reference 

K 

0.042 

This  work 

35-1 

0.064 

This  work 

32.0 

0.061 

Petrucci 

31-5 

0.070 

Petrucci 

28.4 

0.091 

This  work 

24.4 

0.103 

Petrucci 

18.9 

0.095 


Matheson 


25.9 
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ABSTRACT 

A  divalent  cation -selective  electrode  was  used  to  study  the  ion 
association  thermodynamics  of  MgSO.,  and  CaSC>4  at  temperatures  between 
15  and  )5  °  C  and  at  various  ionic  strengths.  The  ionic  strength  was 
adjusted  with  sodium  chloride  because  it  is  the  major  salt  in  most  natural 
waters.  The  thermodynamic  association  constants  for  MgS04  and  CaSC4 
compare  well  with  those  determined  from  conductance  studies 

INTRODUCTION 

The  study  of  ion  association  in  solution  has  been  an  active  area  of 
chemical  research  for  decades.  Almost  every  knovn  physico-chemical  method 
has  been  applied  to  this  problem.  Potent lonetrlc  methods  have  not  been 
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particularly  popular  in  the  past  due  to  the  licited  number  of  ions  for 
which  highly  specific  good  quality  electrodes  were  available.  The  present 
availability  of  the  new  generation  of  membrane  and  solid  state  ion  selective 
electrodes  should  change  that. 

In  this  laboratory  we  are  interested  in  the  chenistry  of  natural  waters 
such  as  seawater***  and  geothermal  brines.  Both  in  the  laboratory  and  the 
field,  new  ways  of  treasuring  important  components  of  such  waters  are  badly 
needed.  The  availability  cf  the  new  divalent  cation  electrodes  encouraged 
us  to  examine  their  applicability  to  such  measurements  on  Mg+^  and  Ca4^  ions. 

As  a  test  of  such  utility  we  decided  to  measure  MgSC>4  and  CaSO&  ion  association 
thermodynamics  over  a  range  of  ionic  strength  and  temperature. 

Several  investigators-*^*^  have  carried  out  limited  range  measurements 
on  these  systems  using  various  techniques.  After  the  measurements  reported 
here  were  completed,  Ainsworth**  published  some  CaS(>4  results  obtained  with 
an  Orion  Model  92-20  calcium  electrode.  His  results  at  25°  agree  in  general 
with  those  reported  here.  However,  his  work  was  restricted  to  one  ionic 
strength  fixed  with  CaCl2» 

EXPERIMENTAL 

All  measurements  were  made  with  an  Orion  Model  801  Digital  pH  Meter 
equipped  with  an  Orion  Model  92-32  Divalent  Cation  electrode  and  an  Orion 
90-01  Single  Junction  reference  electrode.  Measurements  were  made  in  a 
large  volume  jacketed  glass  all  whose  contents  could  be  controlled  to  +_  0.  05°  C. 

Reagent  grade  MgCl2»  MgSO^,  CaCl2»  CaS04  and  NaCl  were  used  to  prepare 
stock  solutions.  All  solutions  were  analyzed  using  standard  ion  exchange 
and  EDTA  titration  techniques. 

In  each  sec  of  measurements  the  first  step  is  the  calibration  of  the 
electrode  pair  using  standard  solutions  of  known  activity.  The  standards 
for  MgSO^  and  CaS(>4  were  MgCl2  and  CaCl2  respectively.  The  activities  of 
these  solutions  are  well  known?*®  These  standard  calibration  curves  gave 
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the  Nernst  slope  of  T9. 5  cv  at  25°.  The  calibration  was  checked  with  standard 
solutions  after  each  *>*»t  of  "unknowns’*  to  test  electrode  drift.  Measurements 
where  such  drift  vas  found  were  discarded. 

After  the  initial  calibration  curye  was  run,  "unknown"  solutions  were 
prepared  by  weight  and  placed  in  the  ceil.  A  constant  potential  reading 
was  nornallv  obtained  in  10-15  minutes  with  the  approximate lv  200  cn^  samples. 
The  pH  of  each  "unknown"  was  checked  to  detect  any  hydrolysis  problems. 


RESULTS  A NO  DISCUSSION 


A.  MgS04 


Two  of  the  authors  have  previously  worked  on  the  >'.g£04  system* 
using  a  spectrophotcmetric  technique  and  potent ioraetrv.  However,  these 
oeasurenents  were  carried  out  at  high  ionic  strengths.  In  t  ‘  is  work  we 
have  used  lower  ionic  strengths  and  have  used  NaCl  as  a  medium.  Although 
not  commonly  used  to  control  ionic  strength,  NaCl  is  the  common  salt  in 
most  natural  waters  and,  so,  has  been  used  here. 

Measurements  were  made  at  four  ionic  strengths  and  three  temperatures 
(15°,  25°,  35°)  in  the  NaCl  media.  In  addition,  MgS04  measurements  were  made 
at  25°  at  seven  concentrations  in  pure  water.  Tables  1,  II  present  sample 
data  at  ionic  strength  0. 108  (25°)  and  0.027  (35°).  Table  III  summarises 
the  pure  water  data.  Table  IV  summarizes  the  NaCl  media  data  at  all  three 
temperatures. 

The  potent lometr ic  readings  are  converted  to  Mg4-*  activity  using 
the  calibration  curve.  The  activities  converted  to  concentration  using  the 
known  activity  coefficient 

*e+J  •  4.U  *2  /  v 

It  Is  then  possible  to  calculate  the  stoichiometric  association  constant 

c 

The  infinite  dilution  associatlcn  constant  is  related  to  the  stoichiometric 
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TABLE  IV 

SIWARY  OF  STOICHIOMETRIC  ASSOCIATION 
CONSTANTS  OF  >:gS04  IN  AQUEOUS  NaCl 


T 

°C 

I 

M 

x± 

I* 

Kc 

LIT1 

15 

0.07470 

0.4432 

0.2732 

32.17 

0.1075 

0.4100 

0.3278 

27.76 

0.1402 

0.3840 

0.3744 

24.42 

0.2462 

0.3285 

0.4962 

16.06 

25 

0.07430 

0.4455 

0.2725 

36.68 

0.1084 

0.4100 

0.3292 

28.92 

0.2404 

0. 3860 

0.3746 

26.18 

0.2478 

0.3275 

0.4978 

16.10 

35 

0.02720 

0.5540 

0.1648 

68.04 

0.07460 

0.4435 

0.2730 

37.74 

0.1085 

0.4080 

0.3293 

28.33 

0.1412 

0.3830 

0.3757 

24.72 

0.2480 

0.3270 

0.4979 

16.97 

constant  by  the  equation 

k.  -  kc 

where  we  heve  assumed  the  activity  coefficient  of  the  ion  pair  to  be  one. 
For  the  calculation  of  K*  ve  have  used  the  Davies  equation.1! 

-  lot\  -  0.5**  A  ■  -  ■*] 

Table  V  summarises  the  results  of  the  Investigation.  The  Ka  in  NaCl  is 
distinctly  lower  than  the  Ka  in  pure  water.  This  is  probably  due  to  minor 
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TAELE  V 

ASSOCIATION  CONSTANTS  AND  THERMODYNAMI C 
PARAMETERS  OF  MgS04 . 


T 

Medium 

A.h° 

£cc 

A  s° 

°C 

LM-1 

KCal/Mole 

KCal/Mole 

KCal/Mole-°K 

IS 

NaCl 

160.83  + 

6.00 

1.15 

-  2.92 

C.0141 

25 

NaCl 

170.46  + 

10.16 

1.15 

-  3.09 

0.0142 

h2o 

204.40  + 

20.00 

35 

NaCl 

182.99  + 

19.67 

1.15 

-  3.20 

0.0141 

amounts  of  NaS04  formation  as  well  as  the  inadequacy  of  the  Davies  equation 
in  the  more  complex  NaCl/Mga04  mixture. 

The  Ka  (H2O)  results  agree  reasonably  well  with  the  results  obtained 
by  conductance  (Ka“187)*2  an<j  other  methods.  13 
B.  CaS04 

Measurements  of  the  CaS04  -  NaCl  system  were  carried  out  at  four 
Ionic  strengths  and  three  temperatures.  Tables  VI,  VII  and  VIII  give 
representative  data.  Each  reported  Kc  is  an  average  of  five  separate 
measurements. 

Calculations  were  carried  out  exactly  the  same  as  in  MgS04.  Table 
IX  summarizes  the  Kc  results  while  Table  X  summarizes  the  and  other 
infinite  dilution  results.  Fig.  1  gives  the  van  Thoff  plot  for  K^. 

The  Ka  value  at  25°  is  196.5  LM“*  in  good  agreement  with  the 
conductance  value*^(200  LM”1)  and  the  solubility  val’ue1^(204. 2  LM”1) .  The 
H°  and  &  S°  are  in  good  agreement  with  those  of  Ainsworth. 


CaS(>4  In  Aqueous  NaCl:15°C  +  0.05;  Y+  -  0.5125;  I  -  0.0402. 
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TABLE  IX 


Sunsaary 

of  Association  Constants 

of  CaSO,  in  Aqueous 

NaCl. 

T 

I 

Y± 

Kc 

°C 

LIT1 

15 

0.0402 

0.2006 

0.5125 

50.950 

0.0599 

0.2447 

0.4705 

41.380 

0.0800 

0.2826 

0.4375 

30.580 

0.0998 

0.3159 

0.4156 

25.956 

25 

0.02060 

0.1435 

0.5860 

69.630 

0.04075 

0.2018 

0.5110 

54.600 

0.06080 

0.2465 

0.4690 

42.380 

0.08010 

0.2830 

0.4360 

37.300 

0.10010 

0.3164 

0.4150 

31.240 

35 

0.02007 

0.1416 

0.5950 

87.000 

0.04000 

0.2001 

0.5140 

64.610 

0.06000 

0.2450 

0.4680 

49.370 

0.08000 

8.2827 

0.4375 

38.680 

0.09990 

0.3160 

0.4155 

32.490 
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TABLE  X 

Association  Constants  and  Thcrcodvnamic  Parameters  of  CaSO,. 

H 
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T 

All0 

Ac° 

A  S° 

°c 

L>rl 

KCal/Mol e 

KCal/Molc 

KC.il/Mole'°K 

15 

172. 74  +  14.18 

2.20 

-  2.964 

+  0.0179 

25 

196.51  +  7.53 

2.20 

-  3.144 

+  0.0179 

35 

221.21  +  18.91 

2.20 

-  3.322 

+  0.0179 

Our  results  show  that  the  specific  ion  electrode  is  very  capable 
of  giving  good  thcnnodynaraic  data  if  proper  caro  is  taken.  The  main 
limitation  is  the  need  to  make  statements  about  activity  coefficients  in 
rather  complex  systems.  This  same  problem  will  also  intrude  in  natural 
water  work. 
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